Field-induced local moments around nonmagnetic impurities in metallic cuprates 
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We consider a defect in a strongly correlated host metal and discuss, within a slave boson mean 
field formalism for the t — t' ~ J model, the formation of an induced paramagnetic moment which 
is extended over nearby sites. We study in particular an impurity in a metallic band, suitable for 
modelling the optimally doped cuprates, in a regime where the impurity moment is paramagnetic. 
The form of the local susceptibility as a function of temperature and doping is found to agree well 
with recent NMR experiments, without including screening processes leading to the Kondo effect. 
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I. INTRODUCTION 

The remarkable character of disorder effects in low- 
dimensional, strongly correlated materials, has been em- 
phasized in recent work on these systems (for a review, 
see Ref. HI). Doping a Mott insulator usually involves 
atomic substitutions which generate random electric po- 
tentials in the material and frequently structural changes 
as well; these defects induce large scale perturbations 
very different from analogous defects in noninteracting 
systems. Early studies of cuprate high temperature su- 
perconductors (HTSC)2'3.4 led to the discovery that non- 
magnetic point defects (typically Zn, Li defects) enhance 
local AF correlations over a wide range of temperatures 
and dopings. substituting Cu atoms in the Cu02 planes. 
Defects produced by electron irradiation also appear to 
produce very similar physical effects as Zn and Li in many 
caseaii. Nuclear magnetic resonance (NMR) spectroscopy 
revealed the main features of this impurity-driven mag- 
netic polarization in the normal (N) and superconducting 
(S) states of underdoped (UD), optimally doped (OP), 
and overdoped (OD) YBCO samples. In the presence of 
a uniform field B, a staggered magnetization (SM) pat- 
tern due to moments on the Cu(2) ions is formed around 
the impurity, with a spatial extent ^, related to the cor- 
relation length of the pure system. This highly corre- 
lated, dynamic entity produces a paramagnetic polariz- 
ability dx which is Curie-like in the UD system, evolving 
to Curie- Weiss-like behavior 6x ^ {T + Q)~^ in the OP 
to OD ranged. Although these were controlled experi- 
ments on systematic impurity substitutions for Cu, it is 
important to realize that similar magnetic phenomena 
are to be expected to occur for intrinsic disorder due to 
the doping process itself, and may dominate some of the 
low-frequency properties of most cuprate samples, par- 
ticularly in the UD regime. 

Because O increases rapidly with doping, and because 
resistivity measurements show that these defects cause 
very rapid scattering at low T, it has sometimes been in- 
terpreted as a Kondo temperature, enhanced in the pres- 
ence of higher carrier densities capable of screening the 



magnetic moment induced by the impurity^. Several ob- 
servations are at odds with a simple Kondo picture: mag- 
netic and transport signatures are quite strong above O 
in the UD regime; ^ is considerably larger than the lat- 
tice spacing; the amplitude of the SM is much bigger 
than what one would expect from the Friedel response to 
a Kondo screened moment. 

The problem of a single nonmagnetic impurity in a 
correlated host material has already received consider- 
able theoretical attention. In the weak-coupling limit, 
several authors^ii^ii^ modelled the problem with a local- 
ized potential added to a Hubbard Hamiltonian treated 
in a Hartree-Fock approximation. They used NMR and 
transport data for the pure system to assign values to 
the parameters of the model, and obtained good agree- 
ment between theory and experiment for impurity Knight 
shifts (in the N and S states) and resistivities (in the N 
state)!" ^^^Yie Hubbard U was tuned to a value quite close 
to a long range AF instability and if the impurity poten- 
tial was chosen to be nonlocal in the N state. The strong 
coupling, [/ — > oo, limit was considered along two main 
lines. One assumes that a magnetic moment has formed 
as a result of the impurity and the Kondo-screening re- 
sponse of the correlated medium is then studie d^^i^^ . The 
other approach models the pseudogap in the UD regime, 
and finds an induced SM (a spinon boundstate) around 
the impurityi^ii^. The problem has also been studied in 
essentially numerical treatment o^^i^^'^^i^^ . 

In this letter, we provide a semiquantitative solution to 
the problem of a single pointlike, nonmagnetic impurity 
in the OP to OD regimes, with neghgible pseudogap, i.e. 
in the metallic N state of a strongly correlated material 
described by the t — t' — J model. Within a mean field 
slave-boson formalism, appropriate to the U ^ oo limit, 
we derive the set of equations describing the paramag- 
netic moments on the planar Cu sites created by a uni- 
form magnetic field. We give an approximate analytical 
solution to these equations which allows us to capture the 
physics at play: a resonant state is formed, producing a 
spatial SM pattern. Its amplitude, which is related to the 
staggered response of the pure system, can be quite large 
but it decays spatially as one moves a few lattice spacings 
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away from the impurity. The local polarizability Sx has 
a Curie- Weiss form with a that depends sensitively on 
doping. The widths of the spinon and holon bands de- 
crease significantly as one moves towards the impurity, 
suggesting a local near-critical region in the vicinity of 
the defect. These results are qualitatively corroborated 
by a fully self-consistent numerical solution of the equa- 
tions. The formation of these moments and their screen- 
ing by the correlated medium stem from the same set of 
carriers, in contrast with the standard Kondo scenario. 



II. HOMOGENEOUS MODEL. 

Our starting point is the t — t' — J model on a square 
lattice, which is commonly considered to capture the low 
energy physics of the Cu02 plane common to all cuprate 
materials. The additional constraint of non double oc- 
cupancy of the sites is handled via the slave boson for- 
malism in which a projected fermion is represented by 
a product of auxiliary ("slave") fields. The Hamiltonian 
of the impurity-free system in the presence of an applied 
magnetic field B then reads 



,19,20 



<i,j>cr <y 



<i,3> 



It describes strongly correlated fermions c. 



{bi are bosons (holons) and //^ are pseudo-fermions with 
spin a (spinous)) on a square lattice, with hopping am- 
plitudes t [t') between nearest- (next-nearest-) neighbor 
sites and nearest-neighbor antiferromagnetic interactions 
J between spins represented by Si — ^fj ■ a ■ fi. The 
fields are subjected to a local constraint J2a "-itr+^l^i ~ 1 
(riio- = fl^fia) which projects out double occupancy from 
the Hilbert space; it is enforced in the functional form of 
Eq. (P) with La grange multipliers \i 

We use a variant of Ubbens and Lee's^i mean field de- 
coupling scheme, appropriate to the gapless spin liquid 
regime when magnetic solutions are included, and intro- 
duce the order parameters: 



< /il/j^ >= X^J < Sf >= m, 

Q^j E^ < bjh N6, 



< b'jb, >-- 



(2) 



where the last expression implements the local constraint 
on the average, and has been given in terms of the av- 
erage hole doping per site 6. The functional form of the 
resulting Lagrangian for the Bose (6) and for the Grass- 
mann (/) fields and the phase diagram determined in this 
approximation are given in Ref. |21I . 



The homogeneous paramagnetic normal state is ob- 
tained with the choices 



1 + 4 Jxo 



iA, = A, (3) 



yielding a uniform stationary Lagrangian Co{x, Q, A) (xo 
is the non interacting Pauli susceptibility for the renor- 
malized spinon band). Correlations affect the effective 
bandwidths of the particles carrying spin and charge. In 
the homogeneous case, these bands reduce to 



— 2ty.b(cos kx + cos ky) + 4|t^ cos cos ky (4) 



with tf = (J/2)x + tQ, t'j = t'Q', tb = 2tx, and t'^ ^ 
2t'x'- Unprimed and primed variables refer to nearest 
and next nearest neighbor amplitudes, respectively. 



III. SEMIANALYTICAL CALCULATION OF 
LOCAL MAGNETIZATION NEAR IMPURITY 

We now assume the presence a single impurity, e.g. a 
zinc atom. Since has a filled shell, there are no 

spinous and no holons. We model this by adding a term 
to Hamiltonian ([T]) that effectively projects out site 0, 



(5) 



with A —>■ oo. In mean field, the charge and spin sec- 
tors can be studied separately, and we denote by Q'' 
{Ga) the Green functions for the holons (spinous) species 
with the A perturbation, when B is present. If an 
impurity-induced magnetic polarization develops in the 
system, with site dependent magnetizations rrii ^ m 
(see Eqs. (|2l3p ). we need to include processes due to 
this magnetic scattering potential. It can be written as 
V = J^a^t^ ^ ■^J2<ij>ai'^3 ~ "^)o'"-i(T if ^c rcplacc B 
by Bm = B — . The full spinon Green's function Go- 
is then formally given by 



Ga = Ga + GaVaGcr, 



(6) 



which gives a self-consistent set of equations for the mag- 
netizations 



J d^fi^)^^G„{i,i;Lj), (7) 



where / is the Fermi function. In the paramagnetic 
regime, this gives us the linear response to the applied 
field B in the form 



lmJd^^{cJ)GOi^,^■,u;) 



,(8) 
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where both Q and G° (the Green's function of 
the defect-free problem) are spin-independent in zero 



field, 



{mj~m)/m and Mi, 



^-^Efc-ilm /dtj/(w)(^C/(i,fc;w)t/(fc,i;a;)^ (fc and j, i 

and j are nearest-neighbors). The stability of a param- 
agnetic solution requires that all the eigenvalues of the 
matrix M = (-My ) be strictly positive. In order to de- 
termine the Si we need to determine . Since Eq. ([5]) 
describes the removal of the site 0, the solution in the 
case of a rigid band would be 



G°(z,0)G°(0,j) 
G0(0,0) 



(9) 



and a similar form for Q''. However, Eq. ^ shows that 
Xijj Qij ^-I'C site dependent, whereas Eq. ^ holds only 
for the homogeneous system. Preserving self-consistent 
determinations of these parameters implies including 
scattering potentials proportional to Xij ~ X ^-nd Qij — Q 
in the Dyson equations for . Enforcing the non dou- 
ble occupancy constraint also requires special care. A 
full solution for these propagators involves a numerical 
calculation (see below). 

Nevertheless, using perturbation theory and controlled 
approximations, we obtain an analytical solution which 
reveals the nature and main features of the induced po- 
larization. To zeroth order, we use Eq. ^ which al- 
lows us to compute the densities of states and to deter- 
mine Xiji Qij ill Eq. On sites close to the impu- 
rity, the potential Eq. ^ pushes states away from the 
edges of the band (Eq. @) and redistributes those inside 
the band. Holons, which sit primarily at the bottom of 
the band are drastically affected, and Qij is strongly su- 
pressed. Xij almost retains its defect-free value, since its 
main contribution comes from spinous at the Fermi level, 
well inside the band. Beyond a characteristic "healing 
length" these parameters recover their unperturbed val- 
ues. We then use these values of Qij and Xij to generate 
the Green functions ,Q to next order in perturbation. 
We do not iterate the process any further, which implies 
that, whithin the healing length, we do not obtain the 
bond order parameters in a self-consistent manner and 
that the non double occupancy constraint is not enforced 
properly. Yet, this truncation, which allows us to han- 
dle analytically tractable expressions, is not too drastic 
a simplification, for two reasons. One is that, for tem- 
peratures comparable to or larger than O, this healing 
length is quite small, as is seen in Fig. [U which shows 
the local spinon bandwidth at site i, tf{i), as a function 
of the distance from the impurity. The second is that 
at all T, the amplitude of induced staggered polarization 
decays very quickly with r^, and we may consider that 
the system settles back into the unperturbed state for r.^ 
larger than ^ (Fig. [J) of order a few lattice spacing. 

Using these approximations, we solve Eq. ([8]), where 
we consider that the only nonzero Si are for sites i sitting 
up to three shells away from the impurity. We noticed 
(see below) that the integral in the expression for Mij 



is proportional to J/tf{i), and the enhancement close to 
the impurity promotes a tendency towards local moment 
formation, i.e. sizable values of the Si. Far from the im- 
purity, this ratio is much smaller and the magnitude of 
the impurity-induced polarization Si goes to zero. Values 
ascribed to the hopping and correlation amplitudes were 
t = 0.45eV, t' = -OAt, J = O.leV, and the field was set 
to B = 7T. The measured values of the Knight shifts for 
the pure system and their 5 and T variations were well re- 
produced if we assign a value 6 = 0.3 to the hole concen- 
tration at optimal doping. Experimentally, optimal dop- 
ing corresponds to d = 0.15 rather than 0.3. A plausible 
reason which explains this difference is that we are using 
a mean field decoupling. Nevertheless, with our choice of 
parameters we get a value of the homogeneous t f Eq. ^ 
extremely close to that determined in the framework of a 
projected Gutzwiller scheme, where the doping is set to 
Qjg 22,23,24,25 ^ As wc pointed out, the large amplitudes 
of the staggered moments near the impurity appear to 
correlate with the ratios J/tf{i). Indeed, the observed 
reduction of i/ compared to the homogeneous case, for 
sites close to the impurity, has two main impacts. One is 
to create an extended effective scattering potential, and 
this enhances the weight of the staggered Fourier compo- 
nent of the local paramagnetic magnetizations m^. The 
other is to increase the magnetic response^^, since ~ in a 
Stoner-like picture ~ a larger value of J/t/(i) brings the 
system locally closer to a magnetic phase. It is notewor- 
thy that in the range of dopings and temperatures that 
we investigated, the smallest eigenvalue of M decreases 
as one decreases i5 and is always more than one order of 
magnitude smaller than the others, which are of order 1. 
Since it is positive, this confirms that the induced mag- 
netization vanishes in zero field. Its smallness indicates 
a resonant state, close to a transition to a bound state, 
but the accuracy of the calculation does not allow one 
to make a stronger statement. It also shows that in the 
absence of the impurity, where C/ = G° , the only solution 
to Eq. ([5]) is Si = 0, for all i. A numerical inspection 
of the sum over k in the expression of Mij reveals that 
the dominant contribution is obtained when k = i, and 
that the integration of this term over uj is proportional 

to 1/tfil). 

In order to give a functional expression for the stag- 
gered polarization, we have sought to fit the solution of 
Eq. ([HI) with a form 



).9(q ■ Yi 



(10) 



for a site at position = (xi,yi) away from the impu- 
rity. The factor ^(q • r^) = 0.5(cos (tt^x^) -I- cos(7rgj/i)) 
allows us to include both commensurate {q — 0) and in- 
commensurate solutions. We found that the best fit to 
the data was obtained for a commensurate modulation, 
when we chose for f{xi,yi) the (square) lattice version 
of the Bessel function Kq, normalized to a nearest neigh- 
bor distance (Fig. This is not a form which emerges 
analytically from the current theory, but rather one mo- 



tivated by rigorous theories for similar problems in IDi. 
Note that according to Reff27l. apart from the under- 
doped regime, m does not vary significantly with T, so 
one may use the above fitting form either for the Sj or 
for the (rrii — m)/ B. 
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FIG. 1: Local spinon bandwidth tf{i) in units of J, as a 
function of the distance from the impurity site in lattice 
constants at different temperatures, for j = 0.22, 5 = 0.3. 

The relative polarization SiiT) is well represented by 
a Curie- Weiss form C/{T -\- 0), as found in experiment. 
The magnitude of si is strongly enhanced compared to 
what we would have found as a result of a standard 
Friedel oscillation (the solution of Eq. ([8]) when one sets 
J = in the definition of the Mij). Let us emphasize 
once again that these features are direct consequences 
of the correlation term and that they are strikingly 
similar with those found in ID for the case of a non- 
magnetic impurity^. Our analytic solution allowed us 
to determine C, 9 and ^ for 5 = 0.28,0.3,0.32 and the 
results are summarized in the plots of Fig. [2l 



IV. FULLY SELF-CONSISTENT EVALUATION 

One uncertainty in the above discussion involves the 
fact that while the system with impurity is electroni- 
cally inhomogeneous, the slave boson constraint has only 
been enforced globally. To check the accuracy of this 
approximation, we perform real space exact diagonaliza- 
tion of Eq. ([T]), plus the impurity potential ([S]), solved 
together with the self-consistently determined local slave 
boson amplitudeai^'^'^Sii^S. The primary effect of the 
constraint, which we now impose locally, appears to be 
to slave the spatial variation of the holon density to that 
of the spinous, and thus eliminate the unphysical free 
bosonic length scale. In fact, the effects of correlations 
are generally mitigated, e.g. the normalized staggered 
magnetization is also reduced relative to Fig. [21 We find 



FIG. 2: Normalized staggered magnetization s(r) = (m(r) — 
m) / m induced by a nonmagnetic impurity in the t — t' — J 
model in the presence of a magnetic field B of 7 Tesla, where 
m is the magnetization of the homogeneous system induced 
by the field, a) normalized magnetization \s\ near impurity at 
T = 25K for J/t = 0.22 and S = 0.3. Solid line: fit to \s{r)\ oc 
s{l)Ko{r/^)/Ko{l/0 for C = 3. b) T-dependence of nearest- 
neighbor normalized magnetization si for same parameters, 
c) Effective moment C and Curie- Weiss temperature Q vs. 
doping S. d) correlation length ^ vs. T extracted from fit 
illustrated in (a). 



that the results of the fully self-consistent evaluation ap- 
pear to agree qualitatively with those of the semianalytic 
approach, but for a smaller, more realistic doping scale. 

Fig. [3] shows the reduction oi tf in the vicinity of the 
impurity. It is qualitatively similar to that found using 
the analytic approach, (Fig. [T]) but we notice differences 
between the two results. In the fully self-consistent calcu- 
lation, both the spinon bandwidth healing length - which 
has a smaller value than that given by the analytic calcu- 
lation - and tf(i) are temperature independent. This is a 
direct consequence of the enforcement of the constraint. 
Holons are slaved to spinous and the spatial variations of 
Qij and Xij with depend on one single characteristic 
(renormalized Fermi) energy. By contrast, in the ana- 
lytic calculation, holons are treated as quasi- free bosons 
and so the spatial variations of Qij depend on fc^T while 
those of Xij are set by the spinon Fermi energy, which is 
proportional to the homogeneous tf. For experimentally 
relevant temperatures, ksT << tf, and so the spatial 
variations of t/(«) track mainly those of the holons. 

In Fig. m we show results of the full evaluation which 



5 



1.0 
0.8 
0.6 
0.4 
0.2 
0.0 



: I 



■ ■■■ ■ Ml I 



o 



- o — Analytic 5=0.3 T=1 OOK 
■ SC 5=0.15 T=100K 
• SC 5=0.15 T=300K 







FIG. 3: Spinon bandwidth in fully self-consistent evaluation 
as a function of distance from the impurity site at filling S — 
0.15 and T = lOOK (filled squares) and 300K (filled circles). 
Bandwidth from semianalytic calculation at 5 = 0.3 and lOOK 
is shown (open circles) for comparison. 



again reproduce the qualitative aspects of experimental 
NMR results on Zn and Li impurities. In Fig. |4^, we 
show the magnetization on the nearest neighbor site in 
an applied 7T field. The low-temperature upturn of this 
magnetization increases in strength as the doping is low- 
ered. It is important to recall that in the mean field 
treatment of the homogeneous system, there is a tran- 
sition to long range antiferromagnetic order as the tem- 
perature and filling are lowered. Thus the enhanced up- 
turns reflect the approach to this mean field transition, 
the best the mean field theory can do to simulate the 
gradual freezing of spin fluctuations in the underdoped 
phases, as documented, e.g. by NMR, /xSR and neutron 
scattering expts.^^'^^'^i^ii^iS&i^i^ai^. 

In each case, the upturn of the (normalized) magne- 
tization on the nearest-neighbor site has been flt to a 
Curie- Weiss form, shown in the figure. The doping de- 
pendence of the prefactor C and the Weiss temperature 
Q are shown in Fig. respectively. Two types of terms 
control the T dependence of si. One is the Friedel-like 
response found in a normal metal, which is quasi- T- 
independent and thus gives a constant si for T » O. 
The other is the large, staggered response caused by a 
local reduction of tf. It gives the main contribution to 
si, at intermediate T (larger than or comparable to Q). 
For the lower dopings, the proximity to a magnetic phase 
affects the small T behavior. These factors modify the 
Curie- Weiss flt, and hence the ((5-dependent) values of C 
and 8. Experimentallyij^, 8 and C are obtained with 
sizable error bars near optimal doping, since 8 varies 
rapidly with S in that range. Despite these limitations, a 
quite reasonable qualitative a gree ment is found between 
our results and those of Ref. 27 over the range of dop- 
ings where our theory applies. Below optimal doping, the 
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FIG. 4: Results for fully self-consistent evaluation of mag- 
netization from slave boson equations, a) Normalized mag- 
netization si on nearest-neighbor site as function of T for 
values of doping S from 0.14 to 0.30. Solid lines show fits to 
si — C/{T + Q). b) Effective moment constant C and Curie- 
Weiss constant O extracted from fits in a) vs. 5. c) Experi- 
mental data from Refs. [l]|27l. showing the magnetization (Sz) 
vs. distance from impurity r (insert) and correlation length 
^ vs. T extracted therefrom; d) Temperature dependence of 
extracted theoretical correlation lengths ^ vs. T. 



current theory is not valid, since the pairing fleld which 
gives rise to the pseudogap in slave boson mean fleld is 
not present. As in the semianalytic calculations, the 8 
scale is somewhat larger than experiment around optimal 
doping; this may be due to the small pseudogap present 
even at optimal doping which is absent from the present 
theory. The pseudogap, as the superconducting gap it- 
self, is knowi*i to promote bound state formation and 
enhance the Curie behavior found in underdoped sam- 
ples. 

In Figs. Hh) and d), we show the spatial extent of the 
magnetic droplet which forms around the impurity com- 
pared with the results of Ref. [23; here too the agreement 
is fairly good. We note furthermore that the length scale 
extracted here is comparable with the antiferromagnetic 
correlation length of the defect-free system'"'-^^, as found 
explicitly in ID spin chains^. Finally, Fig. [5] shows the 
actual distribution of moments rrii on the various sites. 
Far from the impurity, this distribution tends to a flnite 
value, since —>■ toq as r.; —>■ oo. As the temperature is 
lowered further or the coupling J increased, the magne- 
tization oscillations are enhanced further and the values 
on nearby sites of the same sublattice as the impurity ac- 
tually take on negative values (not shown), as observed 
in experiment- The reasonable agreement that is found 
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FIG. 5: Top: magnetization in real space obtained in fully 
self-consistent evaluation for S = 0.15, B — 7T and T — 
lOOK. Bottom: magnetization m cut through impurity site 
along X direction. 

between the results of the semianalytical and numerical 
calculations (see Figs. [T]-|3]) and experiment suggests that 
our approach contains key ingredients required to cap- 
ture the physical mechanism of moment formation and 
screening in correlated systems. 

V. CONCLUSIONS 

We have shown that a simple theory of a nonmagnetic 
impurity in a correlated host described by the slave bo- 
son mean field representation of the t ^ t' — J model 
can explain the basic features of the measured param- 
agnetic response of Zn impurities in YBCO. This theory 
differs from earlier approaches in that it explicitly treats 
the correlations in the strong coupling limit, yet assumes 
a metallic host suitable for discussion of optimal doping. 
The calculated susceptibility is found to be much stronger 
than the weak Friedel-like response expected for a normal 
metal, due not only to the enhanced background density 
of states in the host, but also to local supression of the ef- 



fective fermionic bandwidth around the impurity. Qual- 
itatively, the impurity carves a hole around itself of size 
roughly the pure AF correlation length, and the response 
is therefore somewhat similar to that calculated earlier 
in models of the pseudogap state —"i^; nevertheless, the 
temperature dependence is Curie- Weiss like, rather than 
Curie like, in agreement with experiment. The doping de- 
pendence is also found to be qualitatively in agreement 
with experiment, albeit with a renormalized doping scale. 
By utilizing a fully numerical treatment of the inhomo- 
geneous slave boson problem, we have shown that the 
need for this renormalization arises primarily from an 
overestimation of the local bandwidth suppression due 
to the global enforcement of the slave boson constraint, 
which leads to an unphysical bosonic lengthscale. When 
the constraint is enforced locally, correlations are weaker 
and closer agreement with the realistic doping scale is ob- 
tained. The good qualitative agreement of the results in 
this work with experiment suggest that the screening of 
the moment reflected in the Curie- Weiss form of the sus- 
ceptibility, which is observed to rise steeply as the system 
is doped, need not be due to many-body effects of tradi- 
tional Kondo type. Instead, it arises from the correlation 
"hole" induced around the impurity by the Hubbard in- 
teraction, and can be captured by relatively simple mean 
field theories which ignore the spin-flip scattering which 
usually leads to Kondo physics. The cross section for 
quasiparticles scattered by the magnetic droplet created 
by the impurity is of course different for spins up and 
down; it is furthermore strongly T-dependent due not to 
Kondo screening but to the temperature dependence of 
the paramagnetic moment, as in ID spin chains. 

In principle, the slave boson approach is capable of 
capturing the entire crossover of the induced moment in 
a correlated host problem, from the metallic regime to 
pseudogap regime. To do this within a single formalism 
would be a useful step towards understanding the effects 
of disorder on the cuprate phase diagram, but requires 
the inclusion of pairing effects. Work along these lines is 
in progress. 

Acknowledgments. PJH and WC were partially 
supported by ONR N00014-04-0060, DOE DE-FG02- 
05ER46236 and by visiting scholar grants from C.N.R.S. 
The authors are grateful to H. AUoul and J. Bobroff for 
many enlightening discussions and clarifications of exper- 
iments. 



1 H. AUoul, J. Bobroff, M. Gabay, and P. Hirschfeld, cond- 

mat XXX (2007). 
^ H. AUoul, P. Mendels, H. Casalta, J.-F. Marucco, and 

J. Arabski, Phys. Rev. Lett. 67, 3140 (1991). 
^ A. Mahajan, H. AUoul, G. Collin, and J. Marucco, Phys. 

Rev. Lett. 72, 3100 (1994). 
* J. Bobroff, W. MacFarlane, H.AUoul, P. Mendels, N. Blan- 

chard, G. Collin, and J. Marucco, Phys.Rev.Lett. 83, 4381 



(1999). 

^ N. Bulut, Phys. Rev. B 61, 9051 (2000). 

« Y. Ohashi, J. Phys. Soc. Jpn 70, 2054 (2001). 

^ N. Bulut, Physica C 363, 260 (2001). 

* Y. Ohashi, Phys. Rev. B 66, 054522 (2002). 

^ J. W. Harter, B. M. Andersen, J. Bobroff, M. Gabay, and 

P. J. Hirschfeld, Physical Review B 75, 054520 (2007). 
° H. Kontani and M. Ohno, Phys. Rev. B 74, 014406 (2006). 



7 



11 



G. KhaliuUin and P. Fulde, Phys. Rev. B 52, 9514 (1995). 1-13, 1997 in Cargese, 1998). 



W. Hofstcttcr, R. Bulla, and D. VoUhardt, Phys. Rev. Lett. 

84, 4417 (2000). 
1^ G. KhaliuUin, R. Kilian, S. Krivenko, and P. Fulde, Phys. 

Rev. B 56, 11882 (1997). 
" R. Kilian, S. Krivenko, G. KhaliuUin, and P. Fulde, Phys. 

Rev. B 59, 14432 (1999). 

D. Poilblanc, D. Scalapino, and W. Hanke, Phys. Rev. B 

50, 13020 (1994). 
^® W. Ziegler, D. Poilblanc, R. Preuss, and W. Hanke, Phys. 

Rev. B 53, 8704 (1996). 
" S. Odashima and H. Matsumoto, Phys. Rev. B 56, 126 

(1997). 

^® S. Odashima, H. Matsumoto, and O. Michikami, Physica 
C 336, 287 (2000). 

A. E. Ruckcnstcin, P. J. Hirschfeld, and J. Appel, Phys. 
Rev. B 36, 857 (1987). 
^° G. Baskaran, Z. Zou, and P. W. Anderson, Solid State 
Commun. 36, 853 (1987). 

M. Ubbens and P. Lee, Phys. Rev. B 46, 8434 (1992). 

22 F. Zhang and T. Rico, Phys. Rev. B 37, 3759 (1988). 

23 C. Shih, T. Loo, R. Edcr, C.-Y. Mou, and Y. Chen, Phys. 
Rev. B 92, 227002 (2004). 

2* P. Anderson, P. Lee, M. Randeria, T. Rice, N. Triverdi, 

and F. Zhang (2004). 
25 P. Lee, N. Nagaosa, and X. Wen, Rev. Mod. Phys. 78, 17 

(2006). 

2® M. Gabay, Physica C 235-240, 1337 (1994). 

2^ S. Ouazi, J. Bobroff, H. AUoul, and W. MacFarlane, Phys. 

Rev. B 70, 104515 (2004). 
2* W. Ziegler, H. Endres, and W. Hanke, Phys. Rev. B 58, 

4362 (1998). 

28 J. Han and D.-H. Loo, Phys. Rev. Lett. 85, 1100 (2000). 
3° Z. Wang and P. Loo, Phys. Rev. Lett. 89, 217002 (2002). 
3^ C. Niedermayer, C. Bernhard, T. Blasius, A. Golnik, 

A. Moodcnbaugh, and J. I. Budnick, Phys. Rev. Lett. 80, 

3843 (1998). 

32 Y. Sidis, C. Ulrich, P. Bourgcs, C. Bernhard, C. Nieder- 
mayer, L. P. Rcgnault, N. H. Andersen, and B. Keimer, 
Phys. Rev. Lett. 86, 4100 (2001). 

33 B. Lake, H.M.ROnnow, N. Christensen, G. Aepph, 
K. Lefmann, D. McMorrow, P.Vorderwisch, P. Smeibidl, 
N. Mangkorntong, T. Sasagawa, et al.. Nature 415, 299 
(2002). 

3* C. Panagopoulos, J. L. Tallon, B. D. Rainford, T. Xiang, 
J. R. Cooper, and C. A. Scott, Phys. Rev. B 66, 064501 
(2002). 

35 H. Kimura, M. Kofu, Y. Matsumoto, and K. Hirota2, Phys. 

Rev. Lett. 91, 067002 (2003). 
3« M.-H. Julien, Physica B 329-333, 693 (2003). 
3^^ C. Panagopoulos and V. Dobrosavljevic, Phys. Rev. B 72, 

014536 (2005). 

3® R. Miller, R. Kicfl, J. Brewer, Z. Salman, J. Sonier, 

F. Callaghan, D. Bonn, W. Hardy, and R. Liang, Phys- 
ica B 374-375, 215 (2006). 

3^^ C. Stock, W. J. L. Buyers, Z. Yamani, C. L. Broholm, J.- 
H. Chung, Z. Tun, R. Liang, D. Bonn, W. N. Hardy, and 
R. J. Birgeneau, Physical Review B 73 (2006). 

*° G. Aepph, T. E. Mason, S. M. Hayden, H. A. Mook, and 
J. Kulda, Science 278, 1432 (1997). 

'^^ P. Bourges, in The gap Symmetry and Fluctuations in 
High Temperature Superconductors, edited by J. Bok, 

G. Deutschcr, D. Pavuna, and S. Wolf (cond-mat/9901333. 
Proceedings of NATO ASI summer school held September 



